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Abstract 

An equipment has been developed at the Heinrich-Hertz-Institut in Berlin 
which within the VLF range allows the measurement of the statistical spectral 
amplitude distribution, the ratio of the spectral amplitude at two different 
frequencies and the statistical group time delay of the phase of the spectral 
amplitude, of atmospherics coming from discrete thunderstorms. The group 
time delay is directly related to the phase of the propagation function within 
the atmospheric waveguide and depends on distance, azimuth, frequency, and 
effective height of the waveguide. Its variations with frequency and time have 
been measured for different thunderstorms. They give evidence of the time 
change of the ionospheric reflection characteristics within the VLF range. 
Since the registrations depend on the azimuth of the incoming atmospherics, 
thunderstorms can be located from the group time-delay measurements. 
Moreover, the statistical treatment of the atmospherics allows a quantitative 
determination of the number of lightning flashes within the thunderstorm. 
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1 . Introduction 


Atmospherics are electromagnetic impulses generated by lightning discharges 
and propagating within the waveguide formed by the earth's surface and the 
lower ionosphere. The spectrum of an atmospheric is rather broad with a 
maximum in the VLF range. The impulse becomes dispersed during its 
propagation in the waveguide, Observation of atmospherics provides the 
following information: 

1. The impulse form of the lightning discharges at the origin. 

2. The number of impulses and impulse series cf the lightning 
discharges as dependent on geographical location and time. 

3. The time-, space-, and frequency- dependent transmission 
function of the aforesaid waveguide formed by earth and 
ionosphere. 

Lightning discharges originate according to meteorological conditions. 

About 100 discharges per second are generated on the earth's surface. There 
are some preferred centers cf atmospherics (e.g. , South America and 
Africa) and preferred times of day of occurrence (e.g. , late afternoon in 
middle Europe) (Reference 1). The impulse series is almost arbitrary but 
shows some repetition frequency (Reference 2). The amplitudes of the dis- 
charges are statistically distributed. 

This nature of atmospherics suggests a statistical spectral method for their 
measurement. Fortunately, the sources of atmospherics usually are 
thunderstorms of limited area. Therefore, a method combining direction 
finding with statistics gives the maximum information when applied to the 
three points mentioned above. 



The Heinrich-Hertz-Institut in Berlin-Charlottenburg contains an equipment 
that measures the spectral amplitude, the ratio of the spectral amplitudes, 
and the group time delay of atmospherics at different azimuths. This equip- 
ment is described in section 2. Section 3 gives a brief survey of the measure- 
ment of: the impulse form at the origin, the location and number of lightning 
discharges, and the amplitude of the transmission function of the waveguide. 
Section 4 interprets the group -time -del ay measurements. The group time 
delay is closely related to the phase of the transmission function and therefore 
complements the information available from atmospherics. 

For information about other methods of atmospheric measurement and their 
results see References 3 and 4. 

2. Description of the equipment 

The equipment performs the following functions : direction finding, measure- 
ment of spectral amplitudes and their ratio, and group time delay, as des- 
cribed below in sections 2.1 through 2.5 and illustrated by Figures 1 through 5. 
These equipments are contained in a compact unit and can be operated manually 
or automatically. In manual operation, a polaroid camera is used. Automatic 
operation is briefly described in section 2.5. 

2 . 1 Direction finder 

The direction-finder aerial system consists of a rod antenna and two fixed 
crossed-loop antennas (Figure 1). An atmospheric arriving at the loop 
antennas stimulates the two resonance amplifiers K and tuned at 11 kHz. 
These amplifiers phase-shift the voltage by ±45 degrees; the difference of 
the resulting voltages is obtained by the following stage. The difference- 
voltage amplitude U does not depend on the azimuth. The phase, however, 
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Figure 1 - Direction finder, block diagram. 












is proportional to the azimuth <p of the incoming atmospheric, when compared 
with the output voltage U of the third amplifier K which is connected with 

KZ Z 

the rod antenna. A phase comparison, therefore, yields a voltage U that 
is proportional to the azimuth of the atmospheric. 

2.2 Measurement of spectral amplitudes 

A resonance amplifier tunable between 5 and 50 kHz, connected with a rod 
antenna, measures the vertical component cf the spectral electric field 
strength of an atmospheric. The rectified output voltage U is fed to the 
X-plates of an oscillograph while the Y-plates are fed by the voltage U 

3 . 

of the direction finder described in section 2.1. The second differentiation 
of the output voltage U gives a voltage that so controls the brightness of 
the X-channel that the U-maximum appears on the oscilloscope screen as a 
bright spot. Maximum U is proportional to the spectral amplitude of the 
field strength, and is thus "plotted" against azimuth on the oscilloscope. 

Figure 2 shows a photograph of the oscillograph screen exposed for five 
minutes. Three thunderstorms of different sizes and directions can be 
separated from the underground noise (which consists of atmospherics from 
distant thunderstorms) . 

2.3 Ratio of the spectral amplitudes 

For spectral-amplitude ratio measurements, the narrow-band amplifier 
described in section 2.2 is replaced by two equal amplifiers tuned at different 
frequencies and combined with a voltage-ratio-measuring unit. The output 
voltage is a dc pulse, whose height is a unique function of the ratio of the 
spectral amplitudes at the two frequencies. This unit is again combined with 
the direction-finder unit and a brightness-control unit, (see Figure 3). Each 
thunderstorm produces an accumulation of spots; each spot is due to a single 
lightning discharge. 
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0 2 4 6 8 

SPECTRAL AMPLITUDE VOLTAGE U 

Figure 2 ” Oscillograph of spectral amplitude voltage U versus azimuth at 
9 kHz (5-minute picture). 
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RATIO OF AMPLITUDES 


Figure 3 ~ Oscillograph of spectral amplitude ratio between 6 and 8 kHz versus 
azimuth (5-minute picture). 
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2.4 Spectral group time deh 

Figure 4 is a block diagram of the equipment for the measurement of group 
time delay. Three resonance amplifiers are tuned at the frequencies 
f , f , and f ; f is an arbitrarily chosen frequency; Af is an arbitrarily 
chosen deviation such that 


f = f - t 
1 m 


f =f + Af 

2 m 


These amplifiers measure field strengths whose phases are <f> <£> and d> 

• j. • ■wa 3 i * *2 * 

They are stimulated by an atmospheric. The output voltages of the two first 
stages are mixed. The resulting mixing frequency has the phase . 


In the third stage, the frequency f is multiplied, which gives a phase 2 Om' 
Both resulting phases are compared in a phase-measuring unit. The output 
voltage U is proportional to the phase difference 


t, l + *2- 24 'm = - A ‘ Aai 


( 2 ) 


gr 




. That is, it is proportional to the group time delay (-At ^ anc ^ to *^ e secon< ^ 

■ gr 

derivative of the spectral phase cf the atmospheric with respect to the 


angular frequency oo = 277 f. 


The voltage U is again combined with the voltage Uc of the brightness - 
control unit and wi fV > voltage U_a of the azimuth of the individual 
atmospheric. This leads to a measurement as in Figure 5, where thunder- 
storms of different azimuths and different group time delays can be clearly 
separated. 


2.5 Automatic operation 

In automatic operation, the equipment uses an automatic camera recorder. 

A programmable switching unit manages a regular program of any six kinds 
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Figure 4 - Group time delay measurement, block diagram. 
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Figure 5 ~ Oscillograph of group time delay versus azimuth at 
6kHz; bandwidth Af = 1.4kHz (5-minutepicture). 


of measurements which is repeated every hour. For the past year, regular 
measurements of spectral amplitudes, of the ratio of the spectral amplitudes, 
and of the group time delay have been performed at three widely spaced 
places in Germany. 

A detailed technical description cf this equipment has been given by Heydt 
(Reference 5). 
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3. Measurement of the spectral electric field 
strength of the atmospherics 

The vertical component of the electric field of a single atmospheric at a 
distance p from its origin can be expressed by the Fourier integral 
(Reference 6) 


E(p,t) 


CO 

J b(p,« i) g(o>) 

— CO 


(p, <*>) + 0(^)3 doj 


( 3 ) 


j\|r 

where p = distanced = time, a» = angular frequency, B(p, o>) = b e = the 
transmission function of the waveguide between earth and ionosphere, 

G(o>) = g = the radiative component of the spectral function of the 
atmospheric at its origin evaluated at distance p= 1km, 


The equipment described in section 2.2 measures a value proportional to the 
amplitude of the Fourier component in Equation 3: 

Ab(p, (o) g((t>). 

The number of atmospherics exceeding a certain threshold S (indicated by a 
vertical line in Figure 2) 

-S/Ag b 

N = N 0 J W(g, ca) dg ~ e ° , (4) 

S/Ab 

where the probability function W(g f oj) of the spectral amplitude g is approxi- 
mated by the analytical function 


i ~s/s n 

W(g,<u) = -r e 


and 


g ((o) = J g W(g, w) dg 


is the mean value of g(a>). 
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Figure 6 shows the result cf counting the number cf dots N at different 

thresholds S (Figure 2) for one thunderstorm at three different frequencies. 

Because of the relation Equation 4, the dots can be approximated by straight 

lines (on log paper) which intersect at the point N at S = 0. N is the total 

o o 

number of atmospherics coming from the thunderstorm within the measured 
time interval. 


The slope of the lines equals 

A = Ag Q (w) b(p, &>). 


Measurements of thunderstorms situated at different distances p allow us to 
separate the value g Q (w) from b(p , w) (Reference 6). Thus, the amplitude b 



Figure 6 - Number of atmospherics N exceeding threshold S of the 
receiver for a thunderstorm located at 800 km east from 
the receiving station, Berlin, at three different frequencies. 
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of the transmission function of the waveguide and the mean amplitude g of 

^ o 

the spectral function of the lightning discharges can be derived as functions 
cf distance and frequency. 

The equipment described in section 2.3 measures the ratio 

b(p, ft» 1 )g(o> 1 ) 

b (P> (o 2 )g((o 2 ) 

of each atmospheric. The scattering of the points within a cluster in 
Figure 3 results from the statistical distribution of the values of g for indi- 
vidual atmospherics and from the spaced area of the thunderstorm in which 
lightning discharges are generated at different distances from the receiver. 


Figure 7 shows the ratio of U measured between 6 and 25 kHz to U measured 
at 8 kHz. The well-known maximum of amplitude around 10 kHz is due to the 
overlapping cf the maximum of the first-order mode and the maximum of g 

o 

in this frequency range. We observe the shift of this maximum to lower 
frequencies with increasing distance of the thunderstorm. 


The group time delay unit (section 2.4) measures the phase difference of the 
Fourier component in Equation 3, 


where 


2 

A $ = $(p, (i) ) + $(p, o> ) - 2 $(p, oj ), 

i p m 

2 fi 2 * 2 

A (^represents — Aca 

6 aT 

<£(0, o>) = \E f (o, m) + ib((o) , 


(5) 


and cd a) anc ^ 40 2 arC an § u ^ ar frequencies defined in Equation 1. The 

difference method is necessary in order to eliminate the large phase value 

2 

k ^ in •'3> (k Q = cd/ c is the wave number in vacuum). A $ is a measure of the 
deviation of group velocity from the velocity of light c. 
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0 10 20 


FREQUENCY f (kHz) 

(b) 

Figure 7 - Ratio U (f kHz)/U (8 kHz) (db)vs frequency for two storms ENE 

and WNW from receiving station, Stockert. Vertical lines indicate 
total extension of clusters in Figure 3. Points are centers of 
clusters. 
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The cluster of points in Figure 5 again results from the scattering of the 
phases of the spectral functions of the single atmospherics and from the 
different distances of the atmospherics within a thunderstorm. A detailed 
examination of the group time-delay measurements follows in the next section. 

4. Interpretation of the group time deity measurements 

The name "group time delay" results from the definition of the time t of 

gr 

the propagation of a wave group. The derivative of the spectral phase in 
Equation 3 with respect to the angular frequency should be zero, i.e. , 

— ( <y t + $) = 0 . 

a a) 

From this, it follows that group time 

t = - a <&/$<£>, 
gr 

ooo 

Comparing this with Equation 5 and replacing the differentials by differences 
gives 

A coAt = - A 2 $ (6) 

gr 

in agreement with Equation 2. 

Group time-delay measurements with a reasonable resolution can be carried 
out in the 5- to 10-kHz frequency range. 

From the theory of VLF propagation, it follows that only the first- and 
second- order modes are of importance within this frequency range and at 
distances larger than 500 km. 


13 



For a theoretical interpretation of the phase -difference measurements we 
approximate the transmission function B in Equation 3 by the vertical electric 
field of a normalized vertical electric dipole (Reference 6). If we consider 
only the modes of first and second order, the phase of this transmission 
function is 



(7) 


where 


X +jY = k (1-S ) ~ k C /2, 
n n o n' o m ’ 


p* =p (1 + h/2R), 


D ? -(Y 2 -Y l)p ^ 


F = 


"XT 


D ~1 - j C 2 
n 4 n 


The are the eigenvalues of the modes of order n. The coefficients of 
their real and imaginary parts are approximately (Reference 6) 


X 


n 


, /ov2 2 b 

_ (n - 1/2) tt _ _e 

nl ,2 h 

2k h 
o 


, and 


n 


, , ..2 2 , 
a. (n - 1/2) TT b 
X c 

2 ~g + — , respectively, 

k h ! 

o 


where h is the virtual height of the waveguide, k is the wave number in vacuum, 

o 

and R is the earth's radius. 
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The quantity 


a.- 

1 


0.34 jg 


0'-l f 2/5 


a-1 


( f in kHz ; h i n k m ) 
(a - 4 ; 6- 65) 


is an ionospheric reflection parameter. The factor a, which gives evidence 
of the increase of reflection during the night, is not very well known. It 
can be used as an additional parameter. £ is the daytime height of the 
waveguide. 


The quantity 



3 ' 10-3 f l/2 


( f in kHz ) 


is a measure of the influence of the earth's finite conductivity cr e 

-3 12 

(a =3-10 mhos/m; e =8.85 • 10 sec/m). 

0 O 

.2 

The difference value A (f> of Equation 5 no longer contains the large value 
k jo* (see Equation 7). We now disregard k jd*. With the approximations 
mentioned above, the spectral phase of an atmospheric is 


where 


-5 3/2 

<E>(f, h, p, a) = j/) (f) + x - 5. 1 • 10 f h x 


+ tan 


-r 


x 


sin(8 x) 


( 1 + e ^ X cos(8 x) ] 



(8) 


(p , h, R in km ; f i n kHz) . 
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If we neglect the influence of the second- order mode and of the finite con- 
ductivity of the earth (a fair approximation in this frequency range, at least 
under day time conditions), Equation 8 reduces to 

$ = 0 + X = ij) + 77 2 C p */8h^O). 


From Equations 5 and 6 we now find the frequency-normalized group time 
delay (in j^sec/kHz) as 


At 


gr 


Af 


2ir 


A 2 $ 


c p* 


Aw 


f = f 


16 f 3 h 2 
m 


(- 


Af _\ 2 

^ HI / 


+ 


lo! Ajt 

277 Af 2 


( 9 ) 


m 


(f and Af in kHz; h and p* in km; c in km/sec). 


The first factor on the right side in Equation 9 is the second derivative of 
the phase of the first-order mode: 


277 



16f 3 h 2 


( 10 ) 


This factor which dominates the phase behavior in the frequency range 

between 5 and 1 0 kHz gives evidence of a simple relation connecting At , 

gr 

distance, frequency, and virtual height. The additional small value 
2 2 

A^ip/Af , which is a pure function of f, can be found by an extrapolation of 
the data At gr (P) to the limit p= 1km. 
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GROUP TIME DELAY — ( (isec/kHz) 


Figure 8 shows the theoretical values of the normalized group time delay 
resulting only from the propagation 


277 


A 2 <f> 

.Aw 2 


A 2 J)_ 
Aw 2 . 


versus distance pwith derived from Equation 8 by day (h = 65 km) and 

by night (h = 85 km). The mid frequency is f = 6 kHz, the bandwith is 

Af = 1.4 kHz, and the parameter is CL- 4. By day, the first mode dominates 

almost exclusively at distances greater than 1000 km; at night the influence 

of the second- order mode becomes small at distances greater than about 

2 2 

4000 km. The dashed lines show values of 277 A x/ Act) . 



DISTANCE p (10 3 km) 


Figure 8 - Theoretical values of group time delay due to the propagation 
waveguide versus distance by day and at night. The dashed 
lines are calculated from the first propagation mode. 
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Theory shows that Equation 9 is a fair approximation of the group time delay 
in this frequency range. 

Figure 9 shows the diurnal variation of the normalized group time delay 
(30-minute means) measured from a prolonged thunderstorm at 6800 km 
distance in a WNW direction from the observing station at Stockert near 
Bonn. Germany. 

The thunderstorm disappeared at 1900 UT, May 10. Because of the greater 
attenuation of the electromagnetic waves in the daytime waveguide, the 
nighttime values of the group time delay are much more easily observed 
than the daytime values. 

In Figure 1 0 the nighttime group time delay values at f = 6 kHz and 
4 &f = 1.4 kHz of 20 prolonged thunderstorms observed between April 20 and 
May 10, 1967, and located mainly in a west direction, are plotted against their 
distance p from the receiving station Stockert. The location of the thunder- 
storms was determined from the TV-cloud pictures of the ESSA 5 satellite. 

A linear relationship between and p is suggested by Equation 10, 

From the slope of the curve in Figure 1 0 we can derive, with the help of 
Equation 9, a nighttime value of the virtual height of the waveguide as 

h N = 70 km, 

which is significantly lower than the nighttime values observed at short 
distances and higher frequencies (References 6 and 7). 
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Figure 9 - Measured group time delay versus time for a thunderstorm 
located at 6800 km WNW from the receiving station Stockert. 




s 


and May, 1967 from thunderstorms mainly located west of the 
receiving station Stockert. The firm line is the theoretical 
value of the first mode with h. T = 70 km. 




Figure 11 shows the ratio 


At r (night) 

R ~~ At (day) 
gr 

between the nighttime minimum and daytime maximum value of the group time 
delay of the same 20 thunderstorms that are examined in Figure 10 plotted 
against At^(night)/Af. The arrows indicate upper limit values. With a 
mean value of R = 0.64 from Figure 11 we obtain a daytime value 

h^ = 56 km 

of the virtual height of the waveguide which is again significantly lower than 
the observed values for short propagation paths and at high frequencies 
(References 6 and 7), but fits well with phase measurements for long 
propagation paths (Reference 8). 



Atqr (NIGHT) , v 

GROUP TIME DELAY “ a ( fisec /kHz) 

Af 


Figure 11 - Ratio Atg r (night) / Atg r (day) versus nighttime group time delay 
measured in April and May, 1967, for different thunderstorms 
located mainly west of the receiving station Stockert. 
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Figure 1 2 shows the group time delay of two distant thunderstorms versus 
l/f*, verifying the inverse third-power law of Equation 10. 

In Figure 10 and Figure 12 the straight lines approximating the measured 

3 

values go nearly through the origin at p= 1km or 1/f -* 0. 

2 2 

This means that the spectral function of the atmospherics (S 0/d f 
in Equation 9) value seems to be very small and thus is difficult to determine. 
Our measurements show that Equation 10 describes sufficiently well the 
phase behavior of VLF waves below 10 kHz. Group time delay measurements 
in this frequency range can therefore serve as an excellent tool for the 
location of distant thunderstorms, 

5. Concluding remarks 

This paper gives a short description of an equipment which measures 
spectral amplitude, amplitude ratio, and group time delay of atmospherics 
versus their azimuth in the VLF range, The amplitude measurements are 
interpreted by a statistical method that gives the amplitude of the propagation 
function of the waveguide as a function of frequency and distance and gives 
the mean spectral amplitude of the lightning discharges at their origin. 

By the amplitude -ratio measurements, the frequency and spatial-dependent 
behavior of the field strength within the VLF range can be determined. 

The group time delay is proportional to the second derivative cf the spectral 
phase of the atmospherics with respect to the angular frequency. Its 
theoretical value 

' At gr = A p/ f3 h 2 (11) 

depends on distance p, frequency f, and virtual height h of the waveguide. 

This relation is verified by measurements of distant and prolonged 
thunderstorms. 
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(a) 

FREQUENCY f (kHz) 



(b) 

FREQUENCY f (kHz) 


10 8 7 6 



3 

Figure 12 ~ Group time delay vs 1/f for two storms ENE and WNW from 
receiving station, Stockert. Vertical lines indicate total 
extension of clusters in Figure 5. Points are centers of 
clusters. (Bandwidth Af = 1.4 kHz.) 
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